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I. ABSTRACT 

The relation between the tunnehng current It and the interaction energy/force E in 
an atomic-scale contact is discussed in a frame of a theoretical model established here. 
According to our model, we predict an existence of two characteristic scaling regimes, where 
the tunneling current is either proportional to the interaction energy It ^ E ot to the the 
square of the interaction energy. It ^ E"^. We show that the existence of a given regime is 
basically controlled by two parameters: (i) the energy degeneracy A and (ii) the hopping 
t between electronic levels involved in the interaction process. In addition, we discuss the 
validity of the Bardeen method to determine the tunneling current in the short tip-sample 
distances. 

II. INTRODUCTION 

The characterization and modification of surfaces and nanostructures is one of the most 
pressing challenges of Nanoscience these days. Scanning Probe Microscopy (SPM) has al- 
ready achieved outstanding performance in this field. In particular. Scanning Tunneling 
Microscopy (STM) [T] and Atomic Force Microscopy (AFM) [[2] are two leading experimen- 
tal tools that have been used for this purposes. The first, sensing the tunneling current 
between a probe and a sample, has been already converted to a standard technique widely 
used to image, characterize and modify objects at the atomic scale. However, universal 
application of this method is limited by the requirement of a conductive sample to detect 
the tunneling current. This serious drawback is not an issue anymore for AFM, which senses 
the force acting between the tip and the sample. This is one of the reasons why the latter 
method has been widely adopted in recent years, when the AFM method became widespread 
in different scientific fields (e.g. biology, chemistry and physics). 

Further proliferation of SPM methods is closely tied to a detailed understanding of un- 
dergoing processes during imaging or manipulation. In particular, the tip-sample interaction 
is a key for a better control and interpretation of measurements. Introduction of a so-called 
qPlus sensor [3j allowed simultaneous detection of both tunneling current and forces. This 
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approach opens up an advanced way to characterize objects at the atomic scale combining 
different detection signals at the same time. Moreover, it provides a direct and precise way 
to study the relation between the chemical force and the tunneling current. This ques- 
tion has already received a lot of attention from both theory [131 EHl El] and experiment 
[SHZl IHl [ini [12]. Recently, simultaneous measurements of both atomic force and tunneling 
current with high precision have been reported fTR \T3\ [T3]. Therefore it is evident we have 
come to the stage at which precise simultaneous measurements of the force and the current 
become the standard. This situation calls for deeper theoretical understanding of how the 
tunneling current and force behave in different interaction regimes. In particular, the ques- 
tion how both the tunneling current and the chemical force are correlated with each other 
has received certain attention and it is still not fully resolved [291 El]- 

Simultaneous acquisition of the electron current and the force between nanoscale objects 
would elaborate our knowledge of their physical and material properties and undergoing 
processes in atomic scale. Indeed, direct comparison of these two quantities measured on 
single gold atom wires during a breaking process has already shed more light on the origin of 
conductance quantization in atomic nanowires p] and the relation between their mechanical 
and transport properties. Application of the technique would also bring valuable information 
about the binding mechanism and electron transport through a single molecule anchored to 
metallic electrodes [151 US [IE]- Very recently this kind of experiment has been already 
realized [T7]. Devices integrating electrical and mechanical functionality on the nanoscale, 
so called NEMS (NanoElectroMechanical Systems) systems are another research area (see 
e.g. [19]) where the simultaneous measurement of the current and force could open new 
perspectives. 

Therefore, to understand the the correlation between force and current in atomic scale, 
we provide here a simple model. In this model (i) we discus the implication of the strong 
interaction regime on the electron current through the tip-sample junction and (ii) we explain 
the relation between the chemical force and the tunneling current in the weakly interacting 
regime. In particular, we will show in this paper that the relation between the tunneling 
current and the chemical force can be used as an indicator of the quantum degeneracy 
between the frontier electronic states of the tip and the sample involved in the electron 
transfer. 
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III. MODEL 



The aim of this paper is to get more insight into the relation between the tunnehng 
current and the chemical force arising during the formation of an atomic contact between 
an SPM probe and a surface atom. We will introduce a simple model, describing both the 
electron tunneling and the interaction energy between two bodies representing the SPM tip 
and the sample. We will denote the two quantum-mechanical systems corresponding to the 
two interacting bodies as a, (3. The systems are described by eigenfunctions (pa, (pp and 
eigenenergies ea, satisfying the Schrodinger equation for the respective decoupled system 
a or /3. Next, we will consider only the interaction between the two outermost atoms on the 
tip and the sample. Indeed, this approximation can be justified by the exponential decay 
of wave-functions outside the surface, which limits the interaction between the two bodies 
to only the outermost atomic layers. In particular, recently reported single-atom chemical 
identification is based on strong locality of the chemical force, almost completely driven by 
the apex atom of the tip and the nearest surface atom [3H1 ES] • The same reasoning is also 
plausible for the tunneling current [32] • 

We thus start with two decoupled systems, a (surface) and (3 (tip), which can be described 
as follows: 



where ea, and (pa, 0/3 are eigenenergies and orthonormalized wave functions of independent 
systems a and (3, respectively. 

IV. TUNNELING CURRENT 

In the first part of this section we will focus on the far distance regime, where both the 
tip and the sample are only weakly disturbed by the interacting potential. In this regime, 
the electron transport through the gap is driven by the tunneling process. Afterwards, we 
will discuss the close distance regime, in which strong coupling between the two systems, 
tip and sample, occurrs. In this regime, three main factors influence the electron transport: 
(i) a contraction of the piezo distance between the tip and sample |37l HI]; (ii) a change 




(1) 



and 
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of local density of states due to the formation of the chemical bond between the tip and 
the sample [3U1 HH]; and (iii) a multiple scattering effect of electrons in the gap |44|. In 
general, there is no strict definition of a tip-sample distance demarcating these two regimes. 
Nevertheless, the onset of the short-range chemical force and the deviation of the I-z curve 
from the exponential form are signatures of the strong interaction regime. 



A. Weakly interacting regime 

If the two systems are brought into contact, an overlap between wave- functions 0q,, (pjs 
arises, giving non-zero probability to transfer an electron from one side to the other and vice 
versa. Then the Schrodinger equation of the coupled system has the following form: 

("2^^' + K(r - Ra) + Vp{r- Rp)^ = e,^,. (3) 

For simplicity, we consider the potential as a sum of potentials V = Va + Vp that satisfy 
the condition Va-Vp = This condition minimizes an error introduced by the perturbation 
theory (for detailed discussion see [32]). In the weakly interacting regime, the tunneling 
current between two electrodes is commonly described by the Bardeen approach [20j within 
the second order perturbation theory: 

h = ^^Y.K,?^{eo.-e,-eV\ (4) 



or in an integral form: 



^^ = ^1 K,lV(e-el^)p,(6)d6, (5) 



where the tunneling matrix T^^ is expressed via the Bardeen surface integral in terms of 
unperturbed wave functions of the original systems (pa,4>i3- 

Ta,^ = /^('^^ V*^" - 0.V0pds . (6) 
B. Strongly interacting regime 

When the two systems are brought into a close contact and the short-range chemical 
bond between the tip and the sample is formed, the wave functions of both systems are 
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FIG. 1: Diagrammatic view of the electron propagation through the tunnel gap showing contribu- 
tions corresponding to different orders of the perturbation theory, up to the 3rd order. 

strongly modified. In this regime, a more sophisticated description of the electron transfer 
going beyond perturbation theory is required. Usually, the formulation in terms of Green's 
functions is the method of choice [1], where all multiple scattering events of an electron in 
the contact region are involved (see fig. [T|. The conductance between the tip and sample at 
an energy level e can be described within the Landauer [10] formalism as follows: 



where T is the transmission function giving an electron propagation probability between 
the tip and sample at an energy e and r and r"*" define the transfer matrix from the left 
to the right electrode and vice versa. Usually the transmission function is evaluated using 
the combination of ab initio calculations in a local orbital basis and the Green's function 
methods (see e.g. [H HJl HSj). In a previous works |25l [261 SI] 5 we have shown that the 
conductance can be expressed in terms of an effective coupling between the tip a and the 
sample (5 such as: 






(8) 



We can express the effective coupling t^^ in a diagrammatic form (see fig. ^ 



t 



a. 




1 -1 



(9) 
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and similarly 

tfai^F) = [l - tp,a9Li^F)ta,^9U^F)] tp,a, (10) 

where ta,i3,ta,i3 are probabilities that an electron will jump between tip (a) and sample 
and vice versa (see fig.[l]) and g'^''^ is the retarded/advanced Green's function of the uncoupled 
tip or sample. 

Eq. ([s]) can be further simplified. First, we will consider the electron transport at the 
Fermi level ep as this is close to the conditions of experimental measurements on metal 
surfaces P^. The expressions in eqs. Q and (|T0|) can be further simplified adopting a so- 
called wide band approximation. In this approximation, the real part of the Green's function 
tends to zero and only the imaginary part of the Green's function dominates. Therefore, 
using a relation pa = ^'^Qaa, we can recast the expression of the conductance as follows: 

Se^ T^%,l3pa{eF)pi3{eF) 

Now, the transmission function is a function of three variables: the hopping ta^^ and the 
density of states of the tip Pjs and the sample pa- Next, we will find suitable expressions 
for these three parameters so as to get more insight into the dependence of the transmission 
function T on them. Based on our previous arguments about the exponential decay of the 
wave function and its strong spatial localization, we will restrict ourselves to the electron 
transport between the surface adatom and the tip apex. 

In order to provide the detailed analysis of the transmission function dependence on 
different factors, we need to define the hopping elements between the frontier orbitals of 
the tip and sample. Several authors |32l |33] made an attempt to find an analytical formula 
for the hopping elements. Here, we use the formula with the exponential dependence on 
distance proposed by J. M. Blanco et al. j26] : 

taA^) = ie-v^, (12) 



z 



where z is the distance between the tip and sample, the exponent m depends on orbital 
symmetry associated with the angular momentum quantum number (for a detailed discussion 
see [2n]) and the variable Wo means the apparent height of the barrier. 

We find an expression for the surface adatom density of states pa via coupling of the 
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adatom to the substrate (as shown in the inset of fig. |2]): 



where F describes a couphng of the electronic state a located on the surface adatom to 
electronic states of the substrate represented by a band k (see inset fig. [2]). The coupling 
term F can be written as: 

r = T^tIM^f) ^ ^ETaMep). (14) 



In other words, the variable F defines the energy width of the a state. Finally, using eqs. ( 14 ) , 



(13) and the relation e ~ e^, we come to a simple expression of the adatom density of states 



Pa in this form: 

Paie) ^ (15) 

To express density of state of tip apex, for the sake of simplicity, we will consider a limit 
case where the density of states of the tip apex and surface adatom are similar, i.e. pp ^ Pa- 
This assumption allows us to express the transmission as a function of a single parameter 7. 
Where the parameter 7 is a ratio between the hopping ta^js and the width of the electronic 
state F 

7 = (16) 



Using eqn. 16 the transmission coefficient T can be expressed in terms of 7 as: 



= — (17) 

(1 + 7^) 

Furthermore, to keep our model as instructive as possible we adopt another condition. We 
will consider only one valence state for each system, represented by a single s-like atomic 
orbital state. We also omit any modification of the surface density of states, i.e., we consider 
Pa or pp to be independent of the energy. Hereafter, 7 is only a number depending on the 



tip-sample distance via the hopping formula eq. (12) and the behavior of the transmission 
function T can be easily studied. Figure [2] shows the dependence of the transmission function 
T on the ratio 7 between the hopping t^^^ and the surface band width F. The optimal 
condition for a transmission of an electron is achieved when the ratio 7 tends to one, i.e. 
ta,i3 ~ F. In addition, we can identify two limiting cases when the transmission function 
diminishes: (i) the tunneling regime in which 7 — )■ as ta,i3 — and (ii) the narrow band 
limit when the coupling parameter F is negligible compared to the hopping ta,j3- 
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FIG. 2: Dependence of the transmission function T on the parameter 7, which is the ratio between 
the coupUng parameter (or reciprocal density of states) F and the hopping ta^/j. The inset shows 
a schematic view of a surface adatom coupUng to the surface. The interaction represented by the 
hopping elements Tq ^ occurs via the localized state a. on the adatom and a surface band k. 

On the other hand, the effect of the multiple scattering on the transmission function 
T can be demonstrated by comparing the traditional Bardeen approach, corresponding to 
the 2^^ order perturbation theory, with an exact solution that includes an infinite series of 



scattering events. The numerator in eq. (11) corresponds to the Bardeen 2°*^ order term. 



while the denominator stems from a sum of the infinite series of scattering events; see eqs. ([O]) 



and (10). Immediately, we see that the relation T = 47^ corresponds to eq. (17) up to the 
2nd QTj-^QTj- other words, fig. [s] reveals the difference between the Bardeen theory and the 
exact solution. In far distances, both curves coincide, having an exponential dependence on 
distance. However, at shorter distances, the Bardeen theory follows the exponential form 
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z-piezo 



FIG. 3: The figure shows the dependence of the tunnehng current on the tip-sample distance 
according to the Bardeen and the multiple-scattering theory. The inset shows the same T-z plot 
but in logarithmic scale, revealing the exponential dependence of the transmission function, or 
current, on the tip-sample distance. 

tending unphysically to infinity. In contrast, the exact solution converges to the ballistic 
regime with unity transmission T = 1. 

The discussion presented above is only a very basic picture of the electron transport be- 
tween the tip and sample. One should always keep in mind, that there are other parameters 
affecting the electron transfer through the tunneling barrier. These include: 

1. contraction of the piezo distance between the tip and sample due to atomic relaxation 
induced by a mutual tip-sample interaction [371 Hi] > 

2. change of the local density of states on both the tip apex and the surface atoms due 
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to formation of the chemical bond between them [30l H5] , 

3. modification of the tunnehng barrier (see e.g. [I3l [35l HE]), 

4. modification of the surface dipole jlTl US], to which the origin of the atomic con- 
trast in Kelvin Force Probe Microscopy on semiconductor surfaces has been recently 
attributed [19]. 



V. THE INTERACTION ENERGY 



In this section we describe the interaction energy arising from an interaction of two weakly 
interacting systems. In general, the total energy Hamiltonian can be expressed as follows: 



H = H^ + Hp + 6V. 



(18) 



The Hamiltonian in eq. (18) can be written using a second quantization formalism: 



(19) 



where ea (a similar equation holds for ep): 



2m, 



(20) 



and the (real- valued) coefficients tap = tf^a- 



ap 



2mf 



V' + Va{r- Ra) + Vp{r- Rp) ijpdf. 



(21) 



Using eq. (21) and (20) introduced above, the interaction Hamiltonian in the orthogonal 



basis set ^ can be written in the following matrix form: 



H 




(22) 



haa 

\ ^Pa <^p 

Let us assume that the total energies and wave-functions of the decoupled systems a, (3 are 
known. To keep our model instructive, we will only consider the one-electron contribution 
stemming from the mutual interaction between the two systems. Hereafter, it is convenient 



to express the total Hamiltonian in terms of the original atomic wave-functions 



of the 
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decoupled systems (here superscript denotes that the matrix elements are evaluated using 
atomic orbitals on each subsystem a,/3, which are not yet orthogonal): 



( uo uo 

"'aa "-alB 

\ "'(3a "lip 



Cq, t 'Jii'aa 
''Pa 



el + 5h% 



(23) 



where the terms 5h"^^ and mean a small change of the atomic levels e°, of the 

subsystems a, /3, defined by eqs. ([T|, ([2|, respectively. Matrices denote hopping 

elements between the two subsystems Note that terms and in eq. (23) are 

defined as : 



and the hopping term t° 



2m. 



-V^ + V^{r-R^) + Vp{r-Rp) 



a/3- 



2m, 



+ \4(r - Ra) + Vp{f- Rp) <Ppdr. 



(24) 



(25) 



Comparing eqs. ([T]) and (24) with eq. (23), we obtain directly an expression for 5h°^^: 



6hL 



Vs(f-R 



)aVp 



P)(Pa 



dr. 



(26) 



It is evident that the integral vanishes if the potential Vis^f— Rf^) is short-ranged. Therefore, 
without presence of a long-range potential (e.g. a dipolar coulombic interaction) we can put: 



Sh° = 



e ~ e . 



(27) 



It is worth commenting that the term 5h" represents the first order term in the perturbation 
theory. Non-orthogonality of the atomic wave functions (pa, 0/3 gives rise to a global matrix 
overlap Q = I + S, where 



S 



Saf: 
1^ SiSa 



(28) 



The new (orthonormal) basis ipi is orthogonalized with both systems through the Lowdin 
transformation [2T] : 

^, = 5:670,. (29) 



The operator 6 can be expanded up to the second order in S as 

e-"2^i--s + -s\ 

2 8 



(30) 
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Our goal is to express the orthogonal Hamiltonian H of the whole system in terms of the 
unperturbed atomic hamiltonian H°. The transformed Hamiltonian matrix H in the new 



orthogonal space by the S expansion via eq. (30) can now be expressed as: 



H = e'5if°e-5 ~ fj _ + 3 2X ^„ _ 1 3 2 

V 2 8 / V 2 8 

H" - hsH° + H"S) + lsH"S + 1{H"S^ + S^H°) + 0(5^), (31) 
2 4 8 



where H° is defined by eq. (23). Introducing eqs. ( p3| and (28) into eq. (31) we can obtain, 
with a little algebra, the desired expression for the matrix elements of H. In the next 
sections, we will discuss the intraatomic (diagonal) and interatomic (off-diagonal) terms 
of the Hamiltonian. Finally, we will find the expression for the interaction energy arising 
from the interaction of the two weakly interacting systems. To do this, we will assume the 
following conditions: 

1. the atomic potentials Va and Vp do not overlap i.e. Va ■ Vjb = 0, 

2. three-center terms are negligible, 

3. the diagonal matrix elements h"^^ and /i^^ in the non-orthogonal (j) basis set approxi- 
mately equal the atomic orbital energies e°„, e^^. 

The last condition is valid if the integral Sh"^^ = j 4>a {yj3{f — R/s)^ 0Q,(ir vanishes. In other 
words, if the potential Vg(r — _R^) is short-ranged. 

A. Off-diagonal (interatomic) terms in approximation 



First, we will discus the off-diagonal terms t. From eq. (31), we get 



W - - -^^ajsih^^a + (32) 

If there is a negligible change of atomic levels due to mutual interaction of systems a, /3, e.g. 
represented by long-range coulombic interaction, then using eq. (27) we can write t^p as: 



tap - Kxp - 2'S'a/3(eaa + e^/?)- (33) 
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B. Diagonal (intraatomic) terms in approximation 



The wave-function overlap between the subsystems gives rise to a repulsive term 5^h 
renormalizing the original diagonal elements of each subsystem: 



and 



^ /3 ^ (S 



h° — h° 

"-aa "'Pl: 



h° - h° 

"-PP "-aa 



and using eq. (27) we finally obtain: 



p ^ P 



and 



6% 



PP — -Yl ^aptaP + ^^Sllz 



^P 



p° — P° 

^P 



(34) 



(35) 



(36) 



(37) 



C. The chemical interaction 



So far, the energy interaction terms derived in the previous sections, see eqs. (33) and 



(37), have introduced only the energy corrections due to orthogonalization of wave functions 
of each subsystem. For our discussion about the relation between the tunneling current and 
the chemical force, we also need to take into account the energy contribution emerging from 
the formation of the chemical bond between the frontier states of the tip and sample. We will 
restrict ourselves to the one-electron approximation. Therefore the effective hopping matrix 
element can be expressed via the Bardeen tunneling contribution (for a detailed discussion 
see [22]): 

W = T^p = - <l'c.V(jy})ds. (38) 

Figure |4] shows two limiting cases we will consider for the interaction between the single- 
electron states of the two subsystems: a strongly localized state e/j (tip) interacting with 
a delocalized surface band e^; (b) two degenerate localized states and placed at the 
Fermi level. In the former, the interaction energy due to hybridization of two states is given 
using the standard perturbation theory up to the second order by: 



^^chem ^ y J!^ap^ ^ y _}±ap^ /ggX 
P ^P p ^P 
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FIG. 4: Schematic view of two limiting cases considered in the text. Two interacting atomic 
levels: (a) a strongly localized state (tip) interacting with a delocalized surface band e^; (b) two 
degenerate localized states ea and placed at the Fermi level. 
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5Ef^- = ^iJ^^^iZ^. (40) 

It is important to note that the first order term vanishes here. Namely, in this approximation, 
we do not consider any contribution to the diagonal elements e°^^^ in the perturbation 
potential already expressed in the new orthonormal basis. In particular, this is an important 
point in our discussion of the force vs. current scaling. 

On the other hand, the latter case requires application of degenerate perturbation theory, 
where the lowest order gives: 

5Et^- = Y.\Kp\^Y.Tap, (41) 
5Ef^- = Y.\Kp\^Y.Ta,. (42) 

a a 

An example of this limiting case can be represented by an interaction between two strongly 
localized states, such as two dangling bonds being present on a semiconductor surface as well 
as on the tip. As we will see in the following discussion, it is just the character (localization) 
of the two interacting electronic states which has a fundamental consequence to the relation 
between the tunneling current and the short-range force in atomic contacts. 



VI. RELATION BETWEEN CURRENT & INTERACTION ENERGY 

In this section, we will derive the relation between the tunneling current and the force 
on the onset of the chemical bond formed by the frontier electron states of the tip and 
sample. According to the previous discussion, we already know the tunneling current can 
be expressed using Fermi's golden rule as: 



Itr^{T^f5{e^-ep). (43) 
From this expression it is immediately evident that the tunneling current It in the far distance 



regime is proportional to (T^)^. On the other hand, using eqs. (33), (37), the interaction 
energy is given by: 



al3 



+ 2 ^ ^P^^^f- 



+ E 



chem 



(44) 
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Finally, the one-electron total energy can be expressed as a function of S and T^: 

E^"* ^ S^Ae - ST^ + E'''^"', (45) 
where the term E'^^^'^ is proportional to or T for non- degenerate and degenerate electronic 



states, respectively; see eqs. (41) and (39). 



Three terms appear in eq. (45). The first two terms are a consequence of the orthogo- 
nalization process between 0q,, 0/3 wave functions. The first term, which shifts the diagonal 
elements, decays fast with distance by ~ S"^, therefore we can consider it insignificant at 



far distance. The second term is directly proportional to ST. From eq. (33), we see that 
the overlap S is proportional to the hopping element T; i.e. T S. Using this relation, 
we immediately conclude that the first two terms (i) mostly cancel each other and (ii) are 
proportional to T^. In other words, it implies the interaction energy due to the orthogonal- 
ization process is proportional to the tunneling current as / ~ i?, following the Hofer-Fisher 
prediction [29] : 

E^^ ^ S^Ae - ST"" ~ (T^)'. (46) 



The third term in eq. (45), representing the level hybridization, shows an interesting 
behavior. In more general cases, we have two non-degenerate and/or delocalized states, for 
which the interaction energy is described via the 2"'^ order term of eq. (39): 



EtT ^ (47) 

Here, we immediately find once again that the tunneling current is directly proportional to 
the the interaction energy. 

Nevertheless, the situation changes if two degenerate and localized states are considered. 



Using eq. (41) we have: 



Comparing eqs. (48), (46) and (43), we obtain a non-trivial dependence relating the current 
and the force. Usually, the short-range chemical interaction E'^^^'^ shows a slower increase 
at larger tip-sample distances than the E^"^ energy term. Therefore with a dominance of the 
j^chem f^Qj-j^^ ^]^g relation between the current and the atomic force follows a, I ^ E"^ power 
law as predicted by Chen |311 |3l] . 

Based on our previous analysis, it is evident that the relation between the tunneling 
current and the interaction energy depends strongly on the relative magnitudes of the three 
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a) 




b) 





£p Pp 



FIG. 5: Schematic view of two model systems with 2 (left) and 3 (right) electronic states, respec- 
tively. 



terms in our model described by eq. (45). Indeed, Hofer & Fisher [29] neglected the first 



order term in the interaction energy, justifying this step by the fact that both wave functions 
and potentials of the tip and sample decay exponentially and the perturbation potential of 
tip's wave function is the surface potential and vice versa. This makes the perturbation of 
the diagonal elements negligible. 

Numerical models 



To gain more insight into the relation between the tunneling current and the interaction 
energy as a function of the tip-sample distance, number of states and their energy alignment, 
we will consider a simple numerical model (see Fig. [5]). The model consists of two subsystems 
a, (3 defined by molecular levels with energies and e/3 and coupled by the hopping elements 



rewritten from eq. (12) as: 



tr 



-\z 



(49) 



where m = + + 1, / is the angular momentum quantum number of a given state a or 
(3, z is the distance between the tip and sample and A is the characteristic decay constant 
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related to the effective work function of the system via A ~ The electron occupancy 

of individual levels is defined according to their position with respect to the Fermi level 
(plotted as a dashed line in Fig. |5]). It means (i) the states bellow the Fermi level are 
fully occupied, (ii) the states on the Fermi level are half-filled and (iii) the states above the 
Fermi level are empty. Then the interaction energy is defined as a difference between sum 
of occupied states multiplied by their occupancy of coupled and decoupled systems. 

In the case of a two-level system including only s-like electronic states on both subsystems 
(see fig. [s]), the energy of the coupled system is determined by the secular equation of the 
following matrix: 



y ts~s ^(3 — ^ 

On the other hand, the density of states pa of the a state is given by: 



(50) 



where is the coupling constant to an electrode (see discussion in the section IV B). Then 
the energy of a given electronic level a is written as (see fig. Isl): 



€a = -A/2 - i r„. (52) 

Similarly, the energy of the electronic state is evaluated as: 

= +A/2 - i Tp. (53) 

Here we restrict our discussion to distances at which only a weak chemical force forms 
between the tip and the sample. Therefore, we can neglect a change of the density of states 



given by eq. (51) and the contraction of the tip-sample distance. As we have seen in the 
previous analysis, in this interaction regime the Bardeen approximation is still valid (see 
fig. [3]). Then the conductance at the Fermi energy ep between two levels ea,s and e^j^s can 



be expressed using eq. (43) as: 



Gtiep) = -^Tr{pa{eF)ta,i3Pf3{eF)tf3,a)- (54) 



Solving the secular equation for the two- level system, eq. (50), we obtain the expression for 
the interaction energy (where t means tg-s)'- 



(r, + r^) ^-(r„ + r^)^ + 4(t2 + a + r.r^^ - /mA(r^ - rj) 

-^1,2 — 7. =■= 7. • (55) 
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If the coupling of both levels is equal, i.e. Ta = T^, the interaction energy becomes even 
more simple: 

iV 



Ei,2 = -± VWTA^, (56) 
and immediately, we see that the interaction energy is proportional to t if the condition 



t ^ A is satisfied so that y/t"^ + A"^ = t. By direct comparison with eq. (54) we find the 
relation between the interaction energy E and the current It to be It ~ E'^. In the other 
limiting case, when the relation t < A is satisfied, we use y/f^ + A"^ ~ ^ + ^ ~ 8^ 
obtain E t^. Therefore in this interaction regime, the current and the interaction energy 
(force) are proportional. It ~ E. In conclusion, we have found two regimes according to 
which the interaction energy may be proportional to the current. It ~ E, or have a quadratic 
relation with it. If ~ -E^, in good agreement with our findings in the previous section. 
A more generalized picture, including the dependence on the distance z and the energy 



alignment A, can be obtained by a direct diagonalization of the matrix given by eq. (50). 
The results for different values of the energy alignment A of two s-like states e^^s and e^^^ 
are shown in fig. [61 where the interaction energy and the tunneling current, the latter given 



by eq. (54), are plotted as a function of the distance z. In this particular case, the hopping 
element is written according to eq. (12) as tap = —e~^^, where we used the parameters 



A = 2.828 and Ao = 10 . The density of states is obtained from eq. (51) with the coupling 
constant F = 0.3 for both states a and /3. The logarithmic plots of the interaction energy 
E and the tunneling current It at the Fermi level with different degeneracy A of states ea 
and e/3 reveals the exponential dependence of both variables on the distance. For the sake 
of simplicity we considered the Fermi level ep = and we omitted the constant and 
the voltage (U) dependence when plotting the current It. The results are shown in fig. |6} 
Comparing figs. [6] a) to d) we can clearly distinguish the two scaling regimes It ^ E or 
It ~ E"^, respectively. The existence of a given regime is controlled by the ratio between 
the degeneracy parameter A and the hopping t^/s- Note, in the case of the fully degenerate 
electronic levels (see fig. |6]i)) only the quadratic relation It ~ E"^ occurs. 

In the next step, let's consider a more complex case including 3 levels ea,s, ^a,d and 6/3,^ 
(see fig. |5}d) keeping parameters A, F and Ao the same as in the 2-level model. The model 
mimics surfaces where both s-like and d-like states are present, such as transition and noble 
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2 2.5 3 3.5 4 4.5 5 5.5 

z 



FIG. 6: Logarithmic plots of the tunnehng current It x (y-f^U^ and the interaction energy E 
along the distance for the 2-level system (see fig. [sjon the left). 



metal surfaces. The interaction Hamiltonian has the form: 







t. 



ea^d — e ts-d 

ts-s ts-d ^I3,s ^ ^ ) 



(57) 



where the energy levels are given by eq. (52), (53). The tunneling current is also evaluated 



using eq. (54) as in the case of the 2-level model. 



The dependence of both the current It and the energy E on the distance is shown in 
fig. [rj This model gives the characteristic scaling regimes close to It ^ E or It ~ E'^ again, 



but with more intricate occurrence along the distance z. From eq. (49), it is clear that the 



hopping element tg-s between two s-like orbitals has slower decay than the hopping tg-d with 
distance z. Therefore, in the case of degeneracy between two s-states, the energy interaction 



is dominated by eq. (48) giving rise to the quadratic dependence It ~ E (see fig. Ma,b). On 
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FIG. 7: Logarithmic plots of the tunnehng current It x (^^^U^^ and the interaction energy E 
along the distance for the 3- level system (see fig. [sjon the right). An influence of the degeneracy 
state parameter A on the relation between It and E is seen. The two characteristic scaling regimes 



It ^ E ov It ^ E"^ are found depending on the degeneracy A and the hopping tg^x 
controlled by the distance z. 



s,d) 



the other hand, when the s and d states are degenerate, the system shows more complicated 
behavior, as the scaling regime varies with the distance z (see fig. [7]c,d). In a far distance 
regime, where the relation ts-d ^ -^-^ holds, the interaction energy is mainly driven by 



eq. (48). This leads to the quadratic dependence It ^ E . At closer distances, whith the 



hopping between two s-levels given by ts_s growing faster than the interaction between 
the s-states prevails, tg-d < ts^^, and the linear dependence between the current and force 
is established. It ~ E. Approaching further, the quadratic form It ~ E"^ appears again, 
because the hopping between the s-like levels is much larger than the degeneracy parameter 

s; i.e. tg^g ^g^g. 
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VII. CONCLUSIONS 



In this paper, we have derived a simple theoretical model, which describes the dependence 
of the tunneling current and the interaction energy on distance. Next, we have analyzed 
the relation between the force and current. Two characteristic scaling regimes have been 
identified depending on the degeneracy of frontier electronic states and strength of the 
interaction. Our simple model predicts two different scaling regimes between the tunneling 
current and the interaction energy/force. In particular, if we consider the interaction of two 
non-degenerate and/or delocalized states, the relation / ~ proposed by Hofer & Fischer 
[22] prevails. The existence of this regime has been already proven both experimentally 
P [71 Uni Il3] and theoretically using DFT simulations [131 EH] on metal surfaces. On the 
other hand, in the case of two strongly localized and degenerate levels, such as dangling 
bond states on semiconductor surface, the power law regime J ~ i?^ prevails. Indeed, 
such experimental observation has been already reported on Si surface |5], but the detailed 
understanding of its origin has been missing up to now. Moreover, preliminary results of our 
DFT calculations considering interaction between a simple Si tip and adatom on the Si(lll)- 
(7x7) surface point out the quadratic relation / ~ i?^ on the onset of the chemical force. 
The origin of the quadratic relation / ~ i?^ can be utilized for an advanced characterization 
of the electronic states of the tip and sample based on the electron transport between them. 

In addition, we discussed the tunneling current dependence on distance using the standard 
Bardeen approach and the approach where all multiple scattering events of electrons passing 
through the tip-sample junction are taken into account. At far distances both approaches 
give identical results. However, the Bardeen approach gives unrealistic increase diverging 
to infinity at close distances, while the multiple scattering method converges correctly to 
the ballistic regime. The effect should be taken into account when analysis of STM images 
acquired at low- resistance regime is performed [5U] . 
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